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The role of the Hall effect in the global structure and dynamics 
of planetary magnetospheres: Ganymede as a case study 

J. C. Dorelli,^ Alex Glocer,^ Glyn Gollinson,^ Gabor Toth^ 

Abstract. We present high resolution Hall MHD simulations of Ganymede’s magne¬ 
tosphere demonstrating that Hall electric fields in ion-scale magnetic reconnection lay¬ 
ers have significant global effects not captured in resistive MHD simulations. Consistent 
with local kinetic simulations of magnetic reconnection, our global simulations show the 
development of intense field-aligned currents along the magnetic separatrices. These cur¬ 
rents extend all the way down to the moon’s surface, where they may contribute to Ganymede’s 
aurora. Within the magnetopause and magnetotail current sheets. Hall currents in the 
reconnection plane accelerate ions to the local Alfven speed in the out-of-plane direc¬ 
tion, producing a global system of ion drift belts that circulates Jovian magnetospheric 
plasma throughout Ganymede’s magnetosphere. We discuss some observable consequences 
of these Hall-induced currents and ion drifts: the appearance of a sub-Jovian ’’double 
magnetopause” structure, an Alfvenic ion jet extending across the upstream magnetopause 
and an asymmetric pattern of magnetopause Kelvin-Helmholtz waves. 


1. Introduction 

The rapid reconnection-driven response of magneto- 
spheric convection to changes in the orientation of the ex¬ 
ternal magnetic field is one of the great puzzles of magne¬ 
tospheric physics. The nearly collisionless magnetospheric 
plasma should effectively shield it from changes in the ex¬ 
ternal field. In the resistive mangetohydrodynamics (MHD) 
limit, magnetospheric reconnection should occur in vanish¬ 
ingly thin extended current sheets on time scales much too 
slow [Parker, 1957; Sweet, 1958] to explain magnetic storms 
and substorms. Either there is some anomalous plasma re¬ 
sistivity (e.g., produced by turbulence) enhancing the Sweet- 
Parker reconnection rate, or higher-order terms in the Gen¬ 
eralized Ohm’s Law (e.g., the Hall effect, electron pressure 
anisotropy, or electron inertia) are somehow preventing the 
formation of extended electron-scale current sheets. 

The important role of Hall electric fields in allowing fast 
reconnection to occur in the limit of very small plasma resis¬ 
tivity was only beginning to be appreciated in the late 90’s 
[Biskamp et al, 1995, 1997; Ma and Bhattaeharjee, 1996; 
Shay and Drake, 1998; Shay et al, 1999], culminating in 
the “GEM (Geospace Environment Modeling) Reconnection 
Challenge” results summarized in Birn et al [2001]. While 
there is still debate about the specific physical mechanisms 
involved - e.g., Are the dispersive properties of whistler and 
kinetic Alfven waves playing the dominant role [Mandt et al, 
1994; Rogers et al, 2001]? Are finite Larmor radius effects 
playing a crucial part [Hesse et al, 1999, 2001]? - there is 
general agreement that simply including the Hall term in the 
Generalized Ohm’s Law is a huge improvement over resistive 
MHD, producing reconnection rates that are comparable to 
those observed in full particle-in-cell (PIC) simulations. 

It is less clear how the GEM Reconnection Challenge re¬ 
sults “scale up” to very large systems like planetary magne¬ 
tospheres or stellar coronal active regions, whose character¬ 
istic scales L are many orders of magnitude larger than the 
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Larmor radii or inertial scale lengths over which collision¬ 
less reconnection physics dominates. While there was some 
early evidence suggesting that Hall reconnection produced 
Petschek-like [Petsehek, 1964] configurations that rendered 
the reconnection rate insensitive to the system size [Shay 
et al, 1999], subsequent Hall MHD simulations of magnetic 
island coalescence - in which reconnection is driven by the 
ideal MHD instability of a system of magnetic structures 
with characteristic scale A (the island wavelength) - sug¬ 
gested that the reconnection rate (both the instantaneous 
rate and the coalescence time) should scale like (d^/A)^/^, 
where di is the ion inertial length [Dorelli and Birn, 2003; 
Dorelli, 2003; Knoll and Chaeon, 2006]. Recent electromag¬ 
netic PIC simulations [Karimahadi et al, 2011] of large scale 
island coalescence are consistent with the earlier Hall MHD 
results in some ways (predicting a coalescence time that 
scales like (dj/A)^/^), but they nevertheless support the idea 
that the maximum instantaneous reconnection rate is insen¬ 
sitive to the island size (at least over the limited parameter 
range so far accessible by PIC simulations). 

Moving from two to three dimensions produces a 
much richer reconnection landscape. Breaking the two- 
dimensional symmetry permits plasma instabilities (e.g., 
oblique tearing modes [Galeev et al, 1986; Daughton et al, 
2011]) that would have been suppressed in two-dimensions. 
Recently, Daughton et al [2011] have argued that collision¬ 
less magnetic reconnection in a realistic three-dimensional 
magnetosphere should produce extended electron-scale cur¬ 
rent sheets that in turn become unstable to secondary flux 
rope formation, ultimately resulting in a turbulent recon¬ 
nection layer. 

Unfortunately, high performance computing capabilities 
are not yet advanced enough to permit even Hall MHD sim¬ 
ulations of an Earth-sized magnetosphere; the same disper¬ 
sive waves that may play a role in producing fast reconnec¬ 
tion also make the Courant-Lewy-Eriedrichs time step pro¬ 
hibitively small for explicit codes. Thus, it is not surprising 
that most Earth-scale global magnetosphere simulations to 
date still use resistive MHD, either relying on numerical re¬ 
sistivity or ad hoe current-dependent “anomalous” resistiv¬ 
ity to produce fast reconnection. While it has been known 
since the GEM Reconnection Challenge that one can eas¬ 
ily achieve reconnection rates comparable with observations 
and kinetic simulations by using current-dependent resistiv¬ 
ity models with appropriately dialed free parameters (e.g.. 
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Otto [2001]), it is by no means clear that simply getting 
the right local reconnection rate is the only way that kinetic 
scale physics influences global magnetospheric structure and 
dynamics. In this work, we address this basic question: 
Does the local structure of collisionless magnetic reconnec¬ 
tion (specifically, the Hall current structure within the ion 
inertial region) influence the global structure and dynam¬ 
ics (viz., convection and field-aligned current patterns) of 
planetary magnetospheres? 

While Earth may still be out of reach of present day Hall 
MHD codes - and even Mercury pushes the high perfor¬ 
mance computational envelope - Jupiter’s third Galilean 
moon, Ganymede, provides us with an ideal opportunity 
to begin addressing these questions. The distance between 
the surface of the moon and the upstream magnetopause is 
about a factor of 10 larger than the oxygen ion skin depth in 
the surrounding Jovian magnetosphere (e.g., Kivelson et al 
[2004]), providing us with a dj/L << 1 magnetosphere that 
is still computationally tractable. Ganymede is also unique 
in that the Jovian magnetospheric magnetic field is strongly 
anti-aligned with Ganymede’s dipole, so that magnetopause 
reconnection is strongly driven by nearly steady upstream 
conditions, a situation analogous to the conditions that drive 
enhanced magnetospheric convection and magnetic storms 
in Earth’s magnetosphere. Einally, we have several Galileo 
flybys under essentially steady magnetospheric conditions, 
allowing for a relatively straightforward observational test 
of our simulations. 

Recent high resolution resistive MHD simulations of 
Ganymede’s magnetosphere by Jia et al [2008, 2009] have 
produced a global picture of Ganymede’s magnetic field that 
agrees well with Galileo magnetometer observations for the 
six Ganymede flybys (e.g., reproducing such basic features 
as the location and shape of the magnetopause). While this 
may seem at first glance to support the idea that ion-scale 
physics does not have a significant impact on Ganymede’s 
global magnetic field structure, it is important to note that 
only two of the flybys (G 8 and G28) crossed the upstream 
magnetopause at locations close enough to the reconnection 
site to test predictions about global convection and field- 
aligned current patterns. The G7 and G29 flybys did not 
sample the central tail current sheet (where, as we will see 
below, most of the reconnection-driven convection is con¬ 
fined), and the G 1 flyby was too close to the moon’s surface 
to sample the tail current sheet. Thus, simply comparing 
a simulation to a few Galileo flybys is not sufficient to ad¬ 
dress the impact of ion-scale physics on global magneto¬ 
spheric structure; MHD may get the basic size and shape 
of the magnetopause right and still (as we will demonstrate 
below) fail to correctly capture the global convection and 
field-aligned current patterns. 

In what follows, we describe our use of global Hall MHD 
simulations of Ganymede’s magnetosphere to demonstrate 
that the local ion-scale structures observed in the GEM Re¬ 
connection Ghallenge (and all subsequent kinetic reconnec¬ 
tion simulations) have significant impacts on Ganymede’s 
convection and field-aligned current patterns. Specifically, 
the Hall “out-of-plane” magnetic field quadrupole pattern 
that has become one of the defining characteristics of col¬ 
lisionless magnetic reconnection generates a new system of 
field-aligned currents that map directly from the reconnec¬ 
tion sites all the way down to the moon’s surface, modifying 
the region -1 type current system that supports the Alfven 
wing structure close to the moon. Eurther, the same Hall 
current system produces a new J x B force that acceler¬ 
ates ions to their local Alfven speed out of the reconnection 
plane, introducing large asymmetries to the classic Dungey 
[1961] convection pattern. 

In previous pioneering work, Winglee [1994, 2004] empha¬ 
sized some of the effects described in this paper (in partic¬ 
ular, the impact of Hall physics on the field-aligned current 
systems and the important role of non-MHD perpendicular 


ion drifts). However, the picture that emerged from this 
previous work was incomplete in important ways. The par¬ 
ticle simulations of Winglee [1994] were two-dimensional and 
did not properly capture the impact of non-MHD drifts on 
the three-dimensional convection pattern (as we will demon¬ 
strate below). The Winglee [2004] simulations, while three- 
dimensional, used cell sizes a factor of ~ 8 larger than the 
ion inertial length, implying that ion scale effects were likely 
swamped by numerical dissipation. As we will see below, it 
is crucial that one resolve the ion inertial length with at least 
5 computational cells per ion inertial length in three dimen¬ 
sions to see significant differences between the resistive and 
Hall MHD global convection patterns. 

Previous multi-fluid simulations of Ganymede [Paty and 
Winglee^ 2004, 2006; Paty et a/., 2008] that included the 
Hall effect have achieved resolution comparable to that of 
the simulations described here; however, this previous work 
focused on the effects of ionospheric onflow (separately mod¬ 
eling the hydrogen and oxygen populations) and did not dis¬ 
cuss the impact of Hall-mediated reconnection on the global 
pattern of field-aligned currents and convection. Eurther, it 
is not clear from these simulations how the enhanced iono¬ 
spheric outflow impacted the ability to resolve the ion iner¬ 
tial length: di ~ ^ implying that more oxygen flowing 

out into the magnetosphere makes it more difficult to achieve 
> 5 computational cells per di. In contrast, the present work 
does not explore the role of ionospheric outflow, focusing in¬ 
stead on the global effects of ion-scale reconnection (which 
are much more difficult to resolve in the presence of signif¬ 
icant ion outflow). We plan to explore the effect of iono¬ 
spheric outflow in future work using the multi-fluid models 
developed by Gloeer et al. [2009]. 

2. Simulation setup 

We performed our simulations using the BATS-R-US 
global Hall magnetohydrodynamics code developed at the 
University of Michigan [Powell et al., 1999; Toth et al., 
2008]. The Jovian magnetospheric plasma and field param¬ 
eters were the same as those used by Jia et al. [2008] for 
their simulation of the Galileo G 8 flyby. The Jovian mag¬ 
netospheric plasma and field parameters were the same as 
those used by Jia et al. [2008] for their simulation of the 
Galileo G 8 flyby. Ganymede’s dipole moment was set to 
Mx = -716.8 nT, My = 51.8 nT and Mz = -18.0 nT 
in GphiO coordinates (where X points in the direction of 
the incident Jovian magnetospheric flow, Y points along 
the Ganymede-Jupiter line and Z is parallel to the Jovian 
spin axis). The Jovian magnetic field components were set 
to Bx = OnT, By = —6nT and Bz = —77nT. The 
inflow velocity was Vx = 140 km/sec, Vy = 0 km/sec, 
Vz = 0 km/sec. The Jovian magnetospheric mass density 
and pressure were set to p = 56amu — cm~^ and S.SpPa, 
respectively. 

While BATS-R-US provides the option to run the Hall 
model with several non-ideal contributions included in 
Ohm’s law, E = —V x B/c + J x B/(nec) — Vpe/{ne) + 77 J 
(here, V is the plasma bulk velocity, B is the magnetic field, 
c is the speed of light, J is the current density, n is the 
plasma number density, e is the electron charge, and 77 is the 
plasma resistivity), the results presented here used the much 
simpler Ohm’s law: Eihall = —V x B/c + J x B/( 7 rec). 
Our neglect of the electron pressure gradient and resistiv¬ 
ity terms was motivated by our desire to test, in the global 
magnetospheric context, the two primary results of the GEM 
Reconnection Ghallenge papers (summarized in Birn et al. 
[2001]): 1) EifALL is the minimally complex Ohm’s law re¬ 
quired to produce reconnection rates comparable to those 
observed in full PIG simulations, 2 ) the reconnection rate is 
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Figure 1. Hall effects appear when the dissipation scale is sufficiently smaller than di, so that the 
ion and electron bulk velocities become decoupled in thin di scale current sheets. To achieve this scale 
separation, while still recovering the large scale Alfven wing structure at large distances from the moon, 
we used a nested grid in which the highest resolution of 1/32 Rg was concentrated near the moon (left 
panel). This allowed us to achieve a resolution of 5 — 10 computational cells within the di scale current 
sheet (right panel). Thus, since the numerical scheme concentrates numerical resistivity at the grid scale, 
this strategy achieves the required separation between di and the dissipation scale. 


insensitive to the plasma resistivity when the J x B/ (nec) 
term is included in Ohm’s law. As we will see below, relying 
on numerical resistivity also allows us to more easily achieve 
a significant separation (factor of ~ 5 — 10) between di and 
the dissipation scale. 

The Hall and electron pressure terms introduce whistler 
and kinetic Alfven wave dynamics, respectively, both of 
which alter the structure of the reconnection layer in qual¬ 
itatively similar ways (e.g., Rogers et al. [2001]). Thus, by 
neglecting the pressure gradient term in Ohm’s law, we are 
leaving out potentially important kinetic Alfven wave dy¬ 
namics. The kinetic Alfven wave physics, however, does not 
strongly influence the current sheet structure or reconnec¬ 
tion rate for small guide helds; thus, neglecting it is appro¬ 
priate for the G8 flyby case considered here. Euture simu¬ 
lations of other Galileo flybys where the reconnection guide 
held is moderately large will require inclusion of the elec¬ 
tron pressure gradient term. In any case, simply including 
the Hall term at least captures the de-coupling of electron 
and ion bulk velocities which is completely absent in resis¬ 
tive MHD, and this electron-ion decoupling below the ion 
inertial scale allows magnetic reconnection to occur on an 
Alfvenic time scale that is comparable to that observed in 
resistive MHD with current-dependent anomalous resistiv¬ 
ity. 

Our Hall MHD model also neglects ion pressure 
anisotropies and gyro-viscous effects (e.g., the well-known 
“gyro-viscous cancellation” effect [Hazeltine and Meiss, 
1992] in the momentum equation). While these effects - and 
other hnite Larmor radius (ELR) effects - may also have im¬ 
portant consequences for the structure of Ganymede’s mag¬ 
netosphere, we will see below that the Hall term in Ohm’s 
law already captures globally signihcant ion drifts that are 
absent in MHD . Indeed, we argue below that many of the 
new effects observed in our Ganymede simulations (e.g., the 
asymmetric “double magnetopause” boundary layer) have 
also appeared in multi-fluid simulations [Benna et a/., 2010] 
and hylDrid simulations [Muller et a/., 2012] of Mercury and 
can be explained (at least qualitatively) by two-fluid mag¬ 
netic reconnection physics. 

Eor the parameters described above, the Jovian magne¬ 
tosphere ion inertial length was di ^ 450 km = 0.17 Rg- 
Gonvergence experiments in which we compared results 
at increasingly higher resolutions demonstrated that accu¬ 
rately capturing the Hall effect within the magnetopause 


and magnetotail current sheets required at least 5 compu¬ 
tational cells per di; thus, for the results presented here, 
we chose a computational mesh consisting of several nested 
uniform grids with the innermost grid having a cell size of 
AX = AY = AZ = 1/32 Rg ~ S3 km and dimensions 
Lx = 7Rg,Ly = 4.5Rg,Lz = 4Rg (see Eigure 1). This 
innermost grid was chosen to properly resolve the ion iner¬ 
tial scale in a region containing the upstream magnetopause 
and magnetotail current sheets. The outer grid extended 
to 64Rg to capture the Alfven wing structure far from the 
region of interest in the inner magnetosphere. 

Our strategy was to focus on resolving the di scale, rely¬ 
ing on numerical resistivity r] to break the MHD frozen flux 
constraint. Our reliance on numerical resistivity was made 
necessary by the fact that Hall effects do not become signih¬ 
cant until the dissipation scale Irj is a factor of 5 — 10 smaller 
than di] however, resolving a resistive layer of thickness Irj 
with 5—10 cells when Irj is itself a factor of 5 — 10 smaller 
than di would require 25 — 100 computational cells per di^ 
and this is currently still out of reach with present day com¬ 
puters. Thus, relying on numerical resistivity with oc AX 
allowed us to accurately resolve di scale structure with the 
minimal number of computational cells per di. While this 
approach raises questions about how the results depend on 
the resistivity model, we point out that prodigious numerical 
evidence from Hall MHD simulations strongly suggests that 
the Hall-mediated reconnection rate is insensitive to the dis¬ 
sipation physics (e.g., see the GEM reconnection challenge 
papers [Birn et al, 2001]). Gonsistent with these past re¬ 
sults, we have found in our convergence tests that the di 
scale structures begin to converge when AX < di/5, despite 
the fact that the resistivity is numerical. 

We imposed inner boundary conditions on cells within a 
sphere of radius 1.0 Rg- The mass density and temperature 
at the inner boundary were initially set to 550 amu/cm? and 
20 eU, respectively (consistent with Jia et al. [2008]) and 
then allowed to float (with zero radial derivative) thereafter 
(in contrast to Jia et al. [2008], who hxed these values at the 
inner boundary). While Jia et al. [2008] surrounded their 
inner boundary with a conducting shell (to mimic the iono¬ 
sphere) and imposed zero bulk velocity at the solid surface 
inside this conducting shell (thus allowing magnetospheric 
convection to drive ionospheric convection), we allowed the 
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Figure 2. While the Hall effect does not have much impact on the large scale Alfven wing structure, it 
has a dramatic impact on the pattern of field-aligned currents near the moon. Eigures (a) and (c) show 
the resistive MHD case, while (b) and (d) show the Hall MHD case. Streamlines of the magnetic held 
projected into the plane are shown in blue, and bulk velocity streamlines are shown in green. The black 
dashed lines have slope ± cos“^(Un/VA); where Vin is the ambient Jovian co-rotation velocity, and Va is 
the Alfven speed based on the ambient Jovian density and magnetic held. The lower right panel shows 
the characteristic quadrupolar out-of-plane magnetic held pattern that characterizes Hall reconnection. 
These out-of-plane helds are supported by held-aligned currents that merge into the Alfven wings and 
extend all the way down to the moon’s surface. 


azimuthal components of the bulk velocity to hoat (zero 
radial derivative), while the radial component of the bulk 
velocity was chosen so that how into the moon’s surface was 
absorbed and how out of the moon’s surface was zero. The 
magnetic held at the inner boundary was treated by splitting 
off Earth’s dipole field and solving for the non-potential field 
Bi (e.g., see Tanaka [1994]; Powell et al. [1999]). The total 
field was initialized to Ganymede’s dipole field (i.e., Bi = 0) 
and then Bi was allowed to hoat with zero radial derivative 
thereafter. Given the large uncertainties in the ionospheric 
conductivity - and the complex nature of interaction of pre¬ 
cipitating electrons and ions with Ganymede’s atmosphere 
and icy surface (an interaction that involves both sputter¬ 
ing of neutrals and the production of secondary electrons 
from the surface) - we believe this is a reasonable first step. 


Euture work will incorporate a more realistic model of the 
interaction of magnetospheric plasma with Ganymede’s icy 
surface. 

Before running the Hall MHD case, we ran a resistive 
MHD case with identical boundary conditions but using 
an anomalous resistivity of the following form: n = 770 + 
r]A{J/Jc — 1 ), limited by 0 and r]max = 2 x (here, 

J is the magnitude of the current density; 770 = 1 x 
rjA — 2 X IQpw?/s and Jc = 1 x Am~^ are constants) 
As we will see below, our results for this case are comparable 
to those of Jia et al. [2008] and thus serve as a good starting 
point from which to assess importance of the Hall effect. We 
tested a variety of boundary conditions (not reported here), 
and found little impact on the final results. 

Eor the resistive MHD simulations, the BATS-R-US code 
was run in “local time stepping” mode (in which each grid 
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Figure 3. In resistive MHD (left panel), the magnetospheric convection and field-aligned currents 
(EAC) are symmetric about the Jovian magnetospheric inflow direction, exhibiting the classic Dungey 
convection (green arrows on the density plane and purple arrows on the inner sphere) and region-1 EAC 
(blue-white colorbar on the inner sphere) patterns. In constract, the Hall simulation (right panel) shows 
strong asymmetries in both the convection and EAC patterns. In the Hall simulation, Jovian magneto¬ 
spheric plasma appears to enter Ganymede’s magnetosphere along the sub-Jovian flank, populating the 
plasma sheet with Jovian magnetospheric plasma (in contrast to the resistive MHD case, in which the 
wake region is more or less devoid of Jovian magnetospheric plasma). Note the pronounced bending of 
reconnected field lines (green tubes) out of the reconnection plane in the Hall case. 


cell is advanced at the largest time step possible for numer¬ 
ical stability) until a quasi-steady state was reached; then, 
“time-accurate” time stepping was used. Eor the Hall MHD 
runs, we used two different initial states: 1) local time step¬ 
ping as for the MHD simulations, 2) the final time-accurate 
MHD state. We ran the simulations in time-accurate mode 
long enough to reach a quasi-steady state (~ 10 minutes of 
simulated time was found to be sufficient for Ganymede’s 
small magnetosphere). We obtained similar results in both 
cases, giving us confidence that the final Hall MHD state 
presented below was not influenced by the initial condition. 

3. Results 

Eigure 2 shows an overview of the Hall and resistive 
MHD simulation results. The two upper figures show the 
X (GphiO) component of the bulk velocity for the resis¬ 
tive (left) and Hall (right) runs, showing that the large scale 
Alfven wing structure is not significantly modihed by the 
Hall effect. Reconnection at the upstream magnetopause 
and downstream tail current sheets, however, produces the 
characteristic quadrupolar out-of-plane magnetic field pat¬ 
tern (lower right panel). The out-of-plane field is supported 
by a system of intense field-aligned currents (EAG) that 
merge into the large scale Alfven wings far from the moon, 
but these currents also extend all the way down to the 
moon’s surface. 

Looking at the X — Y plane (Eigure 3), however, reveals 
that the Hall effect introduces significant asymmetries into 
the global convection and EAG patterns. In resistive MHD, 
Jovian magnetospheric plasma flows around the upstream 
magnetopause symmetrically, producing a density cavity de¬ 
void of Jovian magnetospheric plasma in Ganymede’s wake. 
This is the classic Dungey [Dungey, 1961, 1963] convec¬ 
tion cycle, driven by magnetic reconnection at the upstream 
magnetopause and downstream tail current sheets. Asso¬ 
ciated with this symmetric convection pattern is an anti¬ 
symmetric pattern of EAG, most pronounced on the flanks 
and vanishing at the upstream and downstream edges of the 
polar cap. 

In contrast, in the Hall MHD simulation Jovian magne¬ 
tospheric plasma enters Ganymede’s magnetosphere along 
the sub-Jovian flank, filling Ganymede’s wake with Jovian 
magnetospheric plasma. The bending of the reconnected 
field lines (green tubes) out of the reconnection plane - 


which produces the classic quadrupolar out-of-plane mag¬ 
netic field pattern shown in the lower right panel of Eigure 
2 - generates a new system of field-aligned currents that 
are most intense near the upstream and downstream edges 
of the polar cap. The global convection pattern also shows 
large asymmetries, with Jovian magnetospheric plasma flow¬ 
ing into Ganymede’s wake region through the sub-Jovian 
flank. We demonstrate below that all of these effects are 
directly driven by Hall electric fields within the upstream 
magnetopause and downstream tail current sheets and, thus, 
should be present to some degree in all planetary magneto¬ 
spheres (e.g.. Mercury and Earth) for which convection is 
driven by collisionless magnetic reconnection. 

Eigure 4 shows a closer view of the EAG and convection 
patterns in the polar cap. The resistive MHD case shows an 
anti-symmetric pattern of EAG into the moon’s atmosphere, 
with parallel current densities of about 0.05/x — Ajm^ (con¬ 
sistent with those reported by Jia et al. [2008, 2009]). In 
the Hall case, the EAG pattern is distorted by the reconnec¬ 
tion EAG system, and now the most intense current densities 
appear in bands centered at the upstream and downstream 
edges of the polar cap. The polar cap convection pattern 
in the Hall case shows localized acceleration channels near 
the upstream and downstream polar cap boundaries, both of 
which are associated with the bending of the magnetic field 
lines out of the reconnection planes in the magnetopause and 
magnetotail current sheets (labelled “JxB” in the figure). 
We will refer to these acceleration channels as “Harang-like 
discontinuities,” since they are similar in appearance to the 
Harang Discontinuity observed in the evening sector convec¬ 
tion electric field of Earth’s ionosphere. 

The Hall reconnection physics is confined to di scale lay¬ 
ers and, as shown in Eigure 2, has little impact on the large 
scale Alfven wing structure. However, Hall electric fields 
appear to have a dramatic impact on the global convection 
pattern in the X — Y plane. Eigure 5 shows cuts of cur¬ 
rent density magnitude and bulk velocity streamlines (green 
lines and arrows) in the X — Y plane for the resistive MHD 
run (left panel) and the Hall run (right panel). The MHD 
case shows the classic Dungey convection pattern, with tail 
reconnection driving flow jets moonward and tailward from 
an X line oriented more or less along the Y axis about 2Rg 
downstream of the moon. Note the symmetric pattern of 
Kelvin-Helmholtz (KH) waves on the flanks, driven by the 
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Figure 4. The resistive MHD polar cap convection pattern shows the classic symmetric Dungey cycle, 
driven by magnetic reconnection at the upstream magnetopause and tail current sheets. Plasma flows 
uniformly over the polar cap from the upstream side to the downstream side, and the EAC pattern shows 
the familiar relationship between flow shear (responsible for magnetic held line bending in the magne¬ 
tosphere) and current density. In the resistive MHD case, the current densities are most pronounced in 
the hanks. In contrast, the Hall MHD convection pattern exhibits large asymmetries, driven by held 
line bending associated with Hall reconnection. The held line bending produces new J x B forces on the 
ions, generating new “Harang-like” discontinuities in the convection pattern. 


how shear between the reconnection return how and the Jo¬ 
vian co-rotation inhow. 

In the Hall case (right panel of Eigure 5), large scale bulk 
how channels appear, howing in the —Y direction within 
the magnetopause and magnetotail current sheets (magenta 
arrows in the hgure). These jets produce an asymmetric 
pattern of KH waves on the upstream magnetopause, with 
larger amplitude oscillations appearing on the sub-Jovian 
side of the magnetopause. Qualitatively, the asymmetry can 
be explained by the asymmetry in the how shear between 
the incoming Jovian co-rotation how (blue arrows) and the 
upstream magnetopause current sheet jet (magenta arrow). 
Jovian magnetospheric plasma enters Ganymede’s magne¬ 
tosphere through the sub-Jovian hank, returning along a 
thickented sub-Jovian boundary layer to join with the up¬ 
stream magnetopause jet. The resulting high how shear on 
the sub-Jovian side of the magnetopause produces large am¬ 
plitude KH waves there. 

A similar ehect explains the patchy appearance of the 
tail current density in the right panel of Eigure 5: KH waves 
driven by the tail ion jet cause the current sheet to undulate 
in and out of the X — Y plane, producing the appearance of 
current dropouts in the plane. Eigure 6 shows how these KH 
waves develop in the tail current sheet. The arrow glyphs 
show the plasma bulk velocity, painted by magnitude, show¬ 
ing how plasma flowing into the tail current sheet interacts 
with the current sheet jet to produce large flow shear on the 
anti-Jovian side, where the KH wave amplitudes are largest. 


4. Comparison with Galileo magnetometer 
observations 

To test the model predictions described above, we probe 
our simulated magnetosphere with a virtual Galileo probe 
corresponding to the G8 flyby. The G8 flyby crossed the 
upstream magnetopause at about 0.75 Rg above the X — Y 
plane. Since upstream magnetopause reconnection should 
be occurring continuously and more or less steadily (given 
the steady strongly southward JME conditions), the G8 in¬ 
bound magnetopause pass should have been an ideal oppor¬ 
tunity for the Galileo magnetometer (see the overview by 
Kivelson et al. [1998]) to observe the out-of-plane Hall elec¬ 
tric field in the exhaust region above the reconnection site. 
Not surprisingly, previous comparisons with resistive MHD 
by Jia et al. [2008] did not see the expected Hall helds, but 
they did note the presence of large amplitude magnetic field 
fluctuations near the magnetopause crossings. Later, Jia 
et al. [2010] used resistive MHD simulations to interpret the 
observed magnetic field fluctuations as the result of bursty 
reconnection producing large scale flux rope structures, or 
Elux Transfer Events (ETEs), similar to those observed at 
Earth’s magnetopause (see the review by [Elphie^ 1995]). In 
what follows, we argue that these large scale magnetic field 
fluctuations on the inbound crossing are in fact steady state 
structures associated with collisionless reconnection. 

Eigure 7 compares observed (blue lines) and simulated 
(green lines) magnetic field components along the G8 flyby. 
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Figure 5. The Hall effect produces large asymmetries in Ganymede’s convection pattern. The left panel 
shows the classic Dungey convection pattern for the resistive MHD case, in which plasma is accelerated 
both moonward and tailward away from an X line more or less aligned with the Y axis. In the Hall case 
(right panel), fast ion jets within the upstream magnetopause and tail current sheets allow Jovian magne- 
tospheric plasma to circulate throughout Ganymede’s magnetosphere, entering from the sub-Jovian flank 
and returning along a thickened sub-Jovian boundary layer to join with the upstream magnetopause jet. 
This asymmetric circulation pattern produces a corresponding asymmetric pattern of Kelvin-Helmholtz 
waves on the upstream magnetopause. 



Figure 6. The Hall effect produces an ion jet in Ganymede’s magnetotail current sheet. This in turn 
produces a large flow shear on the anti-Jovian side of the current sheet (toward the right in the plot). 
Here, we show a cut plane at X = 2.5 (Y is horizontal and Z is vertical), showing how the resulting 
Kelvin-Helmholtz (KH) waves cause an undulation of the current sheet that appears as small scale cur¬ 
rent density drop-outs in the right panel of Figure 5. The arrows show bulk velocity color coded by 
magnitude. 


The top two panels show the resistive MHD results, and 
the bottom two panels show the Hall results. The left pan¬ 
els show the magnetic field components along the published 
G8 trajectory. While there is good qualitative agreement, 
the simulated probe shows a somewhat larger Bx than ob¬ 
served. This is not surprising, since small differences be¬ 


tween the simulated and actual magnetopause shapes may 
produce observable systematic differences in the magnetic 
field components in the closed field line region (which be¬ 
come distorted near the magnetopause). To explore this 
effect, we introduced a small offset to the published G8 tra¬ 
jectory (+0.05 Rg in the X direction, and +0.05 Rg in the Z 
direction). This offset was chosen to improve the agreement 
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Figure 7. Comparisons of our simulations with the Galileo magnetometer observations for the G8 flyby 
show good overall agreement. The left panels show the comparison using the published G8 flyby coor¬ 
dinates. Green lines show the simulation results, and blue lines show the data. The right panels show 
the results of adding to the simulated probe an offset of +0.05 in the X and Z position components. 
The top panels show the resistive MHD case, and the two bottom panels show the Hall MHD case. 


between simulated and observed Bx components, and the 
offset trajectory will be used in all further discussion below. 

Because of the G8 trajectory crosses the magnetopause 
at a hnite Y location, the reconnection plane is not aligned 
with the GphiO X — Z plane. Thus, the Hall out-of-plane 
held would be expected to produce large perturbations in 
both the X and Y components. However, rotating into a 
“Boundary Normal Coordinate” (BNC) system (where the 
new X unit vector points along the magnetopause outward 
normal, the Z unit vector points northward tangent to the 
magnetopause, and the Y unit vector completes the right- 
handed system) would be expected to eliminate the pertur¬ 
bation in the new X component, leaving only the new Y 
(out of the local reconnection plane) component. This is 
exactly what we observe, as we show now below. 

Figure 8 shows the data transformed into BNC coor¬ 
dinates corresponding to the inbound outer current sheet 
crossing, and Figure 9 shows a three-dimensional view of the 
offset G8 flyby trajectory (green line) through our resistive 
(top panel) and Hall (bottom panel) simulations. We deter¬ 
mined the BNC normal direction by inspection of this plot 
for the inbound magnetopause crossing. The corresponding 
coordinate axes are shown in white. Note that the axis near 
the outbound (upper) crossing is simply translated from the 
inbound (lower) crossing, to simplify the interpretation of 
Figure 8, which plots the components in the single inbound 
current sheet BNC coordinate system. 

The yellow rectangles in Figure 8 show where the simu¬ 
lated probe enters and exits the outer magnetopause current 
sheet, which appears as an obvious reversal in the Bz com¬ 
ponent. The second current sheet crossing is less obvious in 
both the simulations and the data and, therefore, we have 


not attempted to identify it in Figure 8 (though it is obvious 
in Figure 9). For reference, the GphiO magnetic held com¬ 
ponents are reproduced in the top two panels, with resistive 
MHD on the left and Hall MHD on the right. As expected, 
transforming the data into inbound magnetopause BNC co¬ 
ordinates eliminates the negative Bx excursion, leaving only 
an enhanced By excursion. This is consistent with what we 
observe in the Hall simulation (right panels of Figure 8): 
The small negative Bx perturbation in the top right panel 
at about Ib.Shours has been transformed away in the bot¬ 
tom right panel, and the small positive By excursion in the 
top right panel has been enhanced in the lower right panel. 
Similarly, when viewed in the inbound BNC coordinate sys¬ 
tem, a large positive Bx perturbation appears on the out¬ 
bound crossing. We have labelled these perturbations “Hall 
Fields” in Figure 8. 

While these features are somewhat difficult to see in the 
line plots of Figure 8, they appear very clearly in the 3D 
view of Figure 9. The resistive MHD case in the top panel 
of Figure 9 shows that there is no held line bending out of 
the reconnection plane as the simulated probe crosses the 
magnetopause boundary layer. In contrast. Hall case in the 
bottom panel shows clearly that there is signihcant bend¬ 
ing out of the reconnection planes of both the inbound and 
outbound crossings. For the inbound crossing, it is obvi¬ 
ous that the perturbed held is predominantly in the BNC 
Y direciton, while the perturbed held is in the X direction 
for the outbound crossing when viewed in the local BNC 
coordinate system for the inbound crossing. This is exactly 
what is observed in the data. 
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Figure 8. The Hall perturbations become particularly clear in Boundary Normal Coordinates (BNC). 
Here, the left panels show the resistive MHD case, and the right panels show the Hall case, with yel¬ 
low rectangles indicating the boundaries of the outer magnetopause current sheet in the simulations. 
The perturbations marked “Hall Eields” in the lower right panel show the out-of-plane magnetic held 
structure in the reconnection exhausts as they are crossed by the simulation probe. Note that similar 
structures are observed in the actual magnetometer data. 


Another interesting feature of the G8 Galileo magnetome¬ 
ter data is the series of large amplitude waves that appear 
(most clearly in the magnetic held magnitude) as Galileo is 
exiting Ganymede’s magnetosphere on the sub-Jovian side. 
These waves do not appear in the MHD simulations, but 
they do appear in the Hall simulations (compare the up¬ 
per left and upper right magnetic held magnitude traces 
in Eigure 8). Eigure 9 shows that the wave-like feature in 
the simulated outbound crossing is caused by the KH waves 
that are driven by enhanced how shear on the sub-Jovian 
side of the upstream magnetopause (as argued above in the 
context of Eigure 5). To our knowledge, these large ampli¬ 
tude wave-like perturbations observed on the G8 outbound 
magnetopause crossing have not yet been interpreted as KH 
waves; our interpretation is that these KH waves - appear¬ 
ing on the sub-Jovian side but not on the anti-Jovian side 
- are a new feature predicted by our Hall model. They 
are driven by the upstream magnetopause ion jet which, in 
turn, owes its existence to ion-scale structure of the recon¬ 
nection site. Thus, we argue that intermittent upstream 
magnetopause reconnection (and associated hux ropes or 
ETEs) is not required to explain the large amplitude mag¬ 
netic held perturbations observed on the G8 inbound and 
outbound magnetopause crossings. These structures can be 
explained as consequences of steady collisionless reconnec¬ 
tion: 1) the classic out-of-plane magnetic held bending, and 
2) the ion drift that supports the upstream magnetopause 
current sheet. Neither of these effects is present in resistive 
MHD simulations. 


5. Physics of the ion drift belts 

While the appearance of the out-of-plane magnetic held 
is an expected and well-understood consequence of the Hall 
effect, the dramatic distortion of Ganymede’s global convec¬ 
tion pattern was a surprise to us. While the sources of the 
new ion jets are conhned to the upstream and tail current 
sheets, they clearly have a global impact on the structure 
of Ganymede’s magentosphere. The most obvious effect is 
the appearance of a new global system of “ion drift belts” 
that help circulate Jovian magnetospheric plasma througout 
Ganymede’s magnetosphere. The most striking difference 
between the classic Dungey resistive MHD pattern and the 
new drift belt pattern is that Ganymede’s wake should be 
hlled with Jovian magnetospheric plasma in a thin (ion iner¬ 
tial scale) layer around Z = 0. Another important global ef¬ 
fect is the appearance of a signihcantly thickened sub-Jovian 
boundary layer. 

Eigure 10 illustrates our new ion drift belts model of 
global convection in Ganymede’s magnetosphere. The hg- 
ure shows another view of the X — Y cut shown in Eigure 
5, this time annotated with a cartoon of the flow topology. 
The red arrows show the Jovian magnetospheric flow; the or¬ 
ange and magenta arrows show the upstream magnetopause 
and tail current sheet ion jets, respectively; the green ar¬ 
rows show how plasma of Jovian magnetospheric origin cir¬ 
culates through Ganymede’s wake region and then back 
upstream, forming a thickened sub-Jovian “double mag¬ 
netopause” structure (note that Jupiter is somewhere be¬ 
low the bottom of the hgure, co-rotating counter-clockwise). 

























































































X - 10 


DORELLI ET AL.: HALL MHD GANYMEDE SIMULATIONS 



Figure 9. Simulations reveal that Ganymede’s magnetopause has a bifurcated current sheet structure 
at the location of the G8 flyby (green line). In the resistive MHD case, the current sheet is relatively 
smooth, and there are no large out-of-plane magnetic held excursions as the simulated probe crosses the 
boundary layer between the current sheets. In the Hall case, however, the asymmetric pattern of Kelvin- 
Helmholtz (KH) waves is apparent, as are the large out-of-plane excursions of the magnetic held (blue 
held lines). In both panels, the white X — Y axes show the orientation of the inbound magnetopause 
crossing Boundary Normal Goordinate (BNG) system. 


The yellow circles show stagnation points of the how. As a 
parcel of Jovian magnetospheric plasma enters Ganymede’s 
wake region on the sub-Jovian hank (1), it splits, with one 
path merging with the tail ion jet (2, 2’) and the other path 
circulating back upstream (3) where it merges with the up¬ 
stream magnetopause jet and collides with the incoming 
Jovian plasma (a) (generating KH waves in the process). 
Meanwhile, the upstream jet (4) collides with the tail jet 
(2’) at stagnation point (c’), where more KH waves are pro¬ 
duced on the anti-Jovian hank. 

The appearance of the Hall-induced ion drift belts is an 
easily understood consequence of local force balance near the 
reconnection X lines. Eigure 11 illustrates the physics of the 
ion out-of-plane acceleration due to Hall currents in the re¬ 
connection plane. In the left panel, the cyan box is the ion 
diffusion region - the region over which the ion and electron 
bulk velocities become de-coupled. In this region, electrons 
lead the ions into the electron diffusion region (yellow box), 
producing an in-plane current in the positive Y direction 
above the X axis and in the negative Y direction below the 
X axis (here, X and Y now refer to the outhow and inhow 
directions in the reconnection plane, respectively). Simi¬ 
larly, electrons lead the ions out of the electron diffusion 
region, producing an in-plane current in the negative X di¬ 


rection to the right of the Y axis and in the positive X direc¬ 
tion to the left of the Y axis. These in-plane current loops 
support the quadrupolar out-of-plane current system shown 
schematically by the red and blue ovals, and the resulting 
J X B force accelerate ions out of the reconnection plane. 
This Hall acceleration mechanism illustrated schematically 
in the right panel. 

We can estimate the speed of the upstream magnetopause 
ion jet by writing down the Y (GphiO) component of the 
momentum equation: 


pVx 


dVv 

dX 


JxBz 

c 


( 1 ) 


where the Hall current is given by Ampere’s law: Jx ~ 
— {clA'K)dBYldZ^ and By is the amplitude of the out-of- 
plane quadrupole held structure. We approximate equation 
(1) as follows: 


VxVv ^ ByBz 
^ 6i Ai4:7T 


( 2 ) 


where Si and are the thickness and length, respectively, 
of the ion diffusion region. Assuming p is approximately 
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Figure 10. The current sheet ion jets produce a global system of “ion drift belts” that circulate Jovian 
magnetospheric plasma throughout Ganymede’s magnetosphere. The left panel shows the current den¬ 
sity magnitude (colors) and flow streamlines (white lines) in the X — Y plane (in a coordinate system 
for which X points opposite to the Jovian co-rotation flow, Z points along the Jovian rotation axis and 
Y completes the right-handed system). The right panel shows the path of a fluid particle as it circulates 
around the Jovian magnetosphere. After entering Ganymede’s tail along the sub-Jovian flank (1), the 
particle path splits into two paths, one of which (2,2’) merges with the tail ion jet while another (3) cir¬ 
culates back around to the upstream side to merge with the upstream magnetopause jet (4). This return 
path (1-2-3-4) results in the formation of a thickened sub-Jovian shear layer (green arrow) that produces 


a “double magnetopause” structure (dashed black 
with the incoming Jovian magnet ospheric plasma 
Kelvin-Helmholtz waves apparent in the left panel. 

constant, mass conservation implies VxXi ^ VzSi ~ 
where Va = Bz j (since ions leave the ion diffusion 

region at the upstream Alfven speed). Eurther assuming 
that the amplitude of the Hall-induced quadrupole field is 
comparable to the Z component of the field upstream of the 
current sheet (i.e.. By ~ Hz, as observed in the simulations) 
shows that the Hall-induced J x B force accelerates ions in 
the Y direction to approximately the upstream Alvfen speed 
as they cross the magnetopause current sheet: Vy ~ Ua- 

Eigure 12 shows a cut of the Y component of the center- 
of-mass bulk velocity (green curve in the bottom panel) 
across the upstream magnetopause near Y = Z = 0. Using 
the value Bz ~ 50 nT just upstream of the magnetopause 
(see the blue curve in the top panel), and the upstream mass 
density of 56amulcm^, we get an upstream Alfven speed of 
Va ~ 146 km/sec, which compares favorably with the ob¬ 
served peak Vy- Thus, we are confident that the Hall effect 
explains the magnetopause ion jet. Similar analysis also also 
confirms that the tail jet owes its origin to the Hall-induced 
J X B force. 

6. Conclusions and Discussion 

We have used high resolution Hall MHD simulations of 
Ganymede’s magnetosphere to demonstrate that Hall cur¬ 
rents within the magnetopause and magnetotail current 
sheets have a significant impact on the global structure of 
Ganymede’s magnetosphere. Specifically, the Hall effect in¬ 
troduces large asymmetries into the field-aligned current 
(EAG) and convection patterns: 

1. Hall-mediated magnetic reconnection produces a new 
system of field-aligned currents that originate in the recon¬ 
nection layers and extend all the way to the moon’s surface. 


lines). The returning particl path (3) also collides 
at stagnation point (a) to produce the sub-Jovian 


2. Hall-mediated reconnection accelerates ions in the 
magnetopause and magnetotail to the local Alfven speed 
in the Y direction, producing a global system of ion drift 
belts that circulate Jovian magnetospheric plasma through¬ 
out Ganymede’s magnetosphere. 

The appearance of the new system of EAG, with maxi¬ 
mum intensities near the upstream and downstream edges 
of the open-closed boundary, is interesting in light of the 
fact that a similar pattern is observed in Hubble observa¬ 
tions (e.g., Feldman et al. [2000] and McGrath et al. [2013]). 
In contrast, in resistive MHD, the region-1 type EAG actu¬ 
ally vanish at the upstream and downstream edges of the 
polar cap and are enhanced on the flanks. This is diffi¬ 
cult to reconcile with observations if Ganymede’s aurora is 
driven by field-aligned currents (as it is in Earth’s magneto¬ 
sphere). While it is premature to claim that the Hall EAG 
system can explain Ganymede’s aurora, it is encouraging 
that collisionless reconnection can serve as a steady magne¬ 
tospheric generator for EAG in the right locations to explain 
the Hubble emission patterns. Euture work should quantify 
this, but such an effort would require a kinetic model of 
the interaction of the precipitating electrons and ions gen¬ 
erated in the magnetosphere with Ganymede’s atmosphere 
and icy surface (both of which may influence Ganymede’s 
atmosphere). 

Some observable consequences of the ion drift belts in¬ 
clude: a) large out-of-plane magnetic field perturbations as¬ 
sociated with the Hall effect in the upstream magnetopause 
and downstream tail current sheets, b) an Alfvenic ion jet 
extending across the low latitude upstream magnetopause, 
c) a thickened sub-Jovian magnetopause boundary layer 
with distinct outer and inner boundaries, d) an asymmet¬ 
ric pattern of magnetopause Kelvin-Helmholtz waves driven 
by the interaction of the incoming Jovian magnetospheric 
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Figure 11 . The acceleration mechanism responsible for the current sheet ion jets is a simple consequence 
of the Hall helds within the ion diffusion regions. In the left panel, the cyan box is the ion diffusion re¬ 
gion, and the yellow box is the electron diffusion region. The blue and red ovals show the quadrupolar 
out-of-plane magnetic field pattern, supported by the in-plane Hall currents. These in-plane currents 
produce an out-of-plane J x B acceleration, shown schematically in the right panel. 


plasma with the upstream magnetopause ion jet, e) an 
Alfvenic ion jet consisting of Jovian magnetospheric plasma 
extending across the central wake region, f) Harang-like dis¬ 
continuities in the ionospheric convection pattern, and g) the 
presence of Jovian magnetospheric plasma in a thin di scale 
layer (centered on Z = 0 close to the moon in the wake 
region. We have argued in this work that some of these 
consequences (the large out-of-plane magnetic held pertur¬ 
bations and the asymmetric pattern of Kelvin-Helmholtz 
waves) are already apparent in the G8 magnetometer obser¬ 
vations. Some of the other consequences should be testable 
using the other Galileo hybys and including other Galileo 
instruments (plasma and waves), but such an investigation 
is beyond the scope of the present paper. More extensive 




Figure 12. A cut along the X axis crossing the upstream 
magnetopause conhrms that the ions are accelerated to 
the local Alfven speed (in this case, 140 km/sec) in 
the out-of-plane direction. The upper panel shows the 
three components of the magnetic held across the up¬ 
stream magnetopause; the bottom panel shows the three 
components of the bulk velocity. 


comparisons between simulations of other hybys and Galileo 
data will be pursued in future work. 

We note that multi-huid simulations by Benna et al. 
[2010] have produced ion drift belts very similar to those 
found in our Ganymede simulations. However, while Benna 
et aL [2010] invoke electron pressure gradient drifts to ex¬ 
plain the ion belts, our simulations did not include the elec¬ 
tron pressure gradient in Ohm’s law. The similarity of our 
results to those of Benna et al. [2010] suggest that Hall JxB 
forces can explain the ion drift belts. We further specu¬ 
late that the “double magnetopause” structure observed in 
our Ganymede simulations is the same structure observed at 
Mercury by Slavin et al [2008] in the MESSENGER data: 
a thickened dusk-side magnetopause boundary layer driven 
by Hall-mediated reconnection (rather than pressure gradi¬ 
ent drifts, as argued by Muller et al [2012]). Eurther simu¬ 
lation work that isolates the Hall JxB force from pressure 
gradient drifts is required to dehnitively answer this ques¬ 
tion. 

Since the effects observed in our simulations are directly 
driven by magnetic reconnection, they should be present in 
all magnetospheres. The large scale consequences and ob¬ 
servable signatures, however, will depend on a number of 
factors: a) the size of the intrinsic magnetic held (which in- 
huences the size of the magnetosphere), b) the role of mag¬ 
netic reconnection in driving global magnetospheric convec¬ 
tion, c) the nature of the magnetized body’s atmosphere and 
surface (which inhuences how field-aligned currents close). 
We expect the Hall effect to have global infiuence in small 
magnetospheres for which convection is driven primarily 
by magnetic reconnection. Eor example, the Hall effect 
likely plays a negligible role in the large scale structure of 
Jupiter’s co-rotating magnetosphere. Mercury, on the other 
hand, possesses a very small magnetosphere with convection 
driven by magnetic reconnection, so that Hall-mediated re¬ 
connection should produce global field-aligned current and 
ion drift belt structures. Mercury’s surface conductivity is 
very small, however, and it is not clear how the Hall-induced 
field-aligned currents will close. Earth lies somewhere in be¬ 
tween: its magnetosphere is relatively large (several orders 
of magnitude relative to the ion inertial length), and its 
convection pattern transitions from an outer reconnection 
dominated region to an inner co-rotation dominated region. 
Understanding the role of ion-scale reconnection in global 
magnetospheric structure is thus a compelling problem that 
will benefit greatly from future comparative simulation stud¬ 
ies. 
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